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ABSTRACT 

This paper considers the characteristics of magnetic bearings used to support a three mass flexible rotor 
operated at speeds up to 14,000 RPM. The magnetic components of the bearing are of a type reported in the 
literature previously, but the earlier analog controls have been replaced by digital ones. Analog-to-digital and 
digital-to-analog converters and digital control software were installed in an AT&T PC. This PC-based digital 
controller was used to operate one of the magnetic bearings on the test rig. Basic proportional-derivative control 
was applied to the bearings, and the bearing stiffness and damping characteristics were evaluated. Particular 
attention is paid to the frequency dependent behavior of the stiffness and damping properties, and comparisons 
are made between the actual controllers and ideal proportional-derivative control. 

INTRODUCTION 

This paper discusses the development of digital controls for magnetic bearings. An analog controlled mag- 
netic bearing developed by the Rotating Machinery and Controls Laboratory at the University of Virginia has 
been described in several articles [1-4], A PC-based digital controller has been implemented for the same mag- 
netic components and flexible rotor. 

Magnetic bearings are coming into increasing commercial use in rotating machinery. Uses include indus- 
trial applications such as compressors, turbines and pumps, and space related applications such as momentum 
wheels, turbomolecular pumps and attitude control devices. Magnetic suspension devices can also be used in 
other applications such as wind tunnels. A full literature review for all magnetic bearings is beyond the scope of 
this paper - the emphasis here is on digital controls for magnetic bearings. 

Single-axis digital control has been reported by several authors. Hisatani, et al. [5] described the magnetic 
support of a 1 kg vertical rotor supported by a digitally controlled axial active magnetic bearing. They noted that 
the particular digital control algorithm used was successful but not quite as qood as the corresponding analog one. 
Carmichael, et al. [6] reported the use of digital controls in a single axis vertical support for a 1 in diameter ball. 
Scudiere, et al. [7] employed a microprocessor based digital control system to support small steel spheres and 
small vertical rotors with masses up to 90 grams. 
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Digital control of magnetic suspension and balance systems for wind tunnels has been described by 
Britcher, et al. [8,9]. Like magnetic bearings, these suspension systems are inherently unstable. Digital controls 
provide the means to stabilize them with a much greater degree of flexibility than analog controls. The system 
employed for the wind tunnels is much more complicated than the single axis control systems indicated above 
because many axes are involved. A PDP 1 1/84 minicomputer was chosen for implementing the digital controller. 

A theoretical treatment of digital control for radial bearings was given by Schweitzer and Ulbrich [10] but 
not tested. Bleuler and Schweitzer [11] analyzed the use of decentralized control for radial magnetic bearings as 
preparation for the application of microprocessor-based digital controls. Gondhalekar and Holmes [12] developed 
a magnetic bearing for control of transmission shaft vibrations. Traxler, et al. [13] described a radial bearing 
using a microprocessor-based digital control. The microprocessor software consisted of a main program and an 
interrupt handler. Also, a magnetic bearing with digital control has been developed utilizing a rotational interrupt 
[14]. The system was employed on a single-mass rotor supported at one end. Orbit and other vibration plots indi- 
cated successful operation of the bearing. 

The purpose of this paper is to demonstrate the relationships between the analog and digital controllers used 
to stabilize the magnetic bearing and the resulting bearing stiffness and damping properties. Earlier work on this 
test rig [ref. 1] has shown the effect of control algorithms similar to those reported here; this work will extend the 
earlier work, with particular attention paid to the frequency dependence of the magnetic bearing dynamic proper- 
ties. Comparison is made between the actual controllers and equivalent ideal controllers which do not exhibit this 
frequency dependence. Both analog and digital control systems are shown to be effective for controlling vibra- 
tions of the test rotor and for exerting a great deal of control over the dynamic properties of the magnetic bearing. 

NOMENCLATURE 


A/D 

Analog-to-digital convertor 

Kr 

Analog derivative feedback gain constant 


area of air gap perpendicular 

K* 

Position probe sensitivity 

to magnetic flux 

K, 

Loop gain for analog PD controller 

C 

Capacitance 

m 

Mass 

c eq { CO) 

Frequency dependent damping 

N 

Number of coil turns per pole pair 

characteristic of magnetic bearing 

PD 

Proportional -de ri vati ve 

D/A 

Digital-to-analog converter (see ZOH) 

R 

Resistance 

D(z) 

Digital control algorithm 

s 

Laplace variable, s = jc o 

F 

Force 

SR 

Sampling rate 

F ext 

Externally applied force 

T 

Sampling period 

g 

length of the air gap 

W 

Weight 

G 

Nominal undeflected gap length (y=0) 

z 

Z-transform variable, z-e sT 

G(s) 

Continuous domain transfer function 

ZOH 

Zero-order-hold device (see D/A) 

1 bx> I by 

Bias current in magnet coil 

P 

sensitivity of the air gap to shaft 

l px> l py 

Perturbation current in magnet coil 


displacement: 

K a 

Power amplifier transconductance 


dx 

KaiD 

A/D converter gain (counts/volt) 

fk) 

Permeability of free space 

K[)!A 

D/A converter gain (volts/count) 

fir 

Relative permeability 

K d 

First derivative feedback gain constant 


Permeability, | x = po(i r 

K eq m 

Frequency dependent stiffness 

% a 

Power amplifier LPF time constant 

characteristic of magnetic bearing 

X A/D 

A/D conversion time 

K g 

Analog proportional feedback gain 

X comp 

Computation time of digital control algorithm 


constant 

*in 

Input LPF time constant (digital controller) 

K u ,K iy 

Actuator gain of magnetic bearing 
in x and y directions 


Output LPF time constant (different 
for analog and digital controllers) 

K my y > ^ mxx 

Magnetic position stiffness (i, j = x or y) 

x. 

Position probe LPF time constant 

K p 

Proportional feedback gain constant 

<s> 

Magnetic Flux 

Kpert 

Perturbation input gain of power amplifier 

0) 

Frequency (rad/sec) 
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SECTION 1 - Description of Magnetic Bearing 

The magnetic bearings used to support the flexible rotor in this work consist of four electromagnets 
arranged radially around the shaft as shown in Figure 1. Varying the currents in the magnet coils determines the 
forces on the bearing journal (the magnetic laminations), thus determining the static and dynamic properties of 
the magnetic bearing. Maslen, et al. [ref. 15] present the governing equation for a bearing of this type, which for 
small displacements about y = 0, can be linearized as 


where 


and 


Kiy ipy K m yyy 

(1) 

\loA g N 2 I by P 

(2) 

G 2 


PoA g N 2 t f 2 by P 2 
r.-i 

(3) 


K iy and K myy are referred to as the actuator gain and the position stiffness , respectively. 

SECTION 2 - Controller Configuration and Block Diagram 

We now wish to incorporate these force relationships (K myy and K iy ) into a block diagram configuration of 
the magnetic bearing. The equation of motion for the bearing journal is (see Figure 2): 

fext * m y + K m yy}' + Kiyipy (4) 

if we assume that i pyo is chosen so that K iy i m =-W, and that i py has been redefined about the operating point i m 
(allowed here since K iy is a linear function of i py ). 

Taking the Laplace transform gives the steady state solution: 

F ex i — (ms 2 + K my y ) Y (s ) + K ly I py (s) (5) 


Re-arranging this equation gives: 


Y(s) 


1 


ms 2 + K„ 


( 6 ) 


F ext ~~ K,ylpy(s) ,,ij ' ‘'myy 

where F ex , - K iy I py (s) represents the total exogenous force applied to the rotor (see Figure 3). 

Inspection of the characteristic equation of this transfer function reveals that the system is open-loop 
unstable, since K myy < 0. Thus, we must employ active feedback control to stabilize the system. Designating the 
controller as G(s), a controller which relates rotor position y (or x) to actuator current i py (or i px ), we have the con- 
figuration shown in Figure 4. 

The CLTF (closed loop transfer function) for the magnetic bearing system with this feedback is 

Y(s) L 


F ext ms + K myy + K iy G(s) 


K myy < 0 


(7) 


which can also be written: 

Y (s)[ms 2 + K myy + K iy G (*)] = (8) 

We are interested in the relationships between the controller G(s) and the stiffness and damping properties 
of the magnetic bearing. For the single degree-of-freedom (SDOF) bearing model shown in Figure 5, the equa- 
tion of motion is: 

Y(s)[ms 2 + C b s + K b ] = F exl (9) 

Thus, we wish to relate the SDOF stiffness and damping coefficients K b and C b to the feedback controller G(s). 
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This relationship must be made with care, however. It is important to realize that, because the forces generated 
by a magnetic bearing are a function of, among other things, a frequency dependent transfer function, G(s), the 
stiffness and damping properties experienced by the bearing journal (and hence the rotor) will be different for dif- 
ferent radial excitation frequencies. We introduce frequency dependent stiffness and damping coefficients, 
K eq ((o) and C eq {L o), to replace terms K b and C b for realizable controllers. It is important to note that (0 does not 
refer to the rotational frequency of the rotor, but rather the radial vibration frequency of the bearing journal. 

Now, we combine equation (8) and equation (9) as follows: 

Y(s)[ms 2 +K myy + K iy G{s)} = Y(s)[ms 2 + C e ,(CO)j + AT e ,(a>)] (10) 

Cancelling like terms, substituting s =j(0 and separating real and imaginary parts of the equations gives: 

K tq (to) = K myy + K iy Re{G (J(o)j (11) 

C eq ((M) = K iy lm{ - G ^ - (12) 

Thus, for the SDOF model of the magnetic bearing/rotor system, we have a relationship between the bearing 
stiffness and damping characteristics and the controller transfer function. 


SECTION 3 - Analog implementation of PD control 

The Routh criterion [ref. 16] dictates that the system described by equation (7) is open-loop unstable. Sta- 
bilization requires that both proportional and derivative terms be present in G(s), the controller transfer function. 
Without considering stability, earlier works have demonstrated the utility of derivative feedback for reducing 
vibration in rotating systems [ref. 17]. This suggests a transfer function of the form: 

G(s) = K p + K d s (13) 


Combining with the expressions for K eq ( to) and C eq ( to) (equations 1 1 and 12) leads to: 

K f q (co) = K myy + K ly K p 

C eq ((0) = K iy K d 


(14) 

(15) 


In other words, the bearing stiffness is a linear function of K p (the proportional feedback coefficient), and the 
bearing damping is a linear function of K d (the derivative feedback coefficient). 

As it turns out, this controller cannot be implemented perfectly. In particular, an infinite bandwidth dif- 
ferentiator will exacerbate the effects of any extraneous noise in the system; high frequency poles must be 
included in the actual design to avoid this problem. Also, the proximity probes and power amplifiers employed in 
a real control system will have finite bandwidth and will introduce additional poles into the controller transfer 
function, G(s). The analog controller reported in this work is represented by: 

0.01^/fJRiC(0.1/if f )5 + AT.(1 +/?Cs)(l +RiC,j)] 

° (lo) 

K) (1+v)(1+RCj)(1+/?iC,j)(1+W)(1+V) 


Figure 6 shows a block diagram of this controller. Note that for the actual analog controller, the proportional 
feedback gain coefficient is denoted by K g and the derivative feedback gain coefficient is written K r . Values for 
the various system parameters described in this paper are presented in Appendix A. 


SECTION 4 - Digital Implementation of PD Control 

The digital implementation of proportional-derivative feedback is quite similar to the analog one. The prox- 
imity probes and power amplifier used with the analog controller are retained in the digital design, after all, the 
interface between any control system and a physical device must ultimately be analog in nature. The difference 
between the analog and digital systems lies mainly with the core elements of the system: while the analog con- 
troller employs resistors and capacitors to realize continuous differential equations, the transfer function of the 
digital controller results from repeated evaluation of a difference equation. The periodic rate at which the com- 
puter evaluates this difference equation is called the sampling rate. 
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Figure 7 is a block diagram of the digital controller. As mentioned earlier, it shares many features with its 
analog counterpart - namely, the proximity probe interface and the power amplifier. There are, however, notable 
differences. Before a signal can be processed by a difference equation, it must be converted to a numerical value, 
thus requiring the presence of an Analog-to-Digital converter between the position signal input and the digital 
computer. Likewise, the periodic stream of numbers generated by successive solutions of the difference equa- 
tions must be converted back to analog voltage signals by a Digital-to-Analog converter prior to power amplifica- 
tion. 

Proportional-derivative control is achieved digitally by converting the analog PD algorithm [equation (13)] 
to its discrete equivalent. The approximation used here is the backwards difference rule for approximating the 
Laplace variable s. The resulting equivalence is 


where T is the sampling period (inverse of sampling rate) and z = e sT . Substituting for s in equation (13) and tak- 
ing the inverse z-transform yields the difference equation: 


K d/ 

~ K„r n + (r B t „_i ) 


(18) 


where r n and c B are points in the input and output data streams of the digital computer, respectively. Note that the 
backwards difference rule approximates the derivative of the incoming signal by determining the slope of a 
straight line intersecting the most recent two data points in the input stream. Detailed discussion of z-transform 
theory and other approximation rules can be found in Franklin & Powell [ref. 18]. 

The digital controller discussed in this paper was implemented on an AT&T PC6300 personal computer 
(IBM PC/XT compatible). Even though this computer has an 8 MHz system clock and is outfitted with the Intel 
8087 math co-processor, its computational abilities are somewhat limited. Thus, in order to achieve the highest 
possible sampling rate for the controller, the software that actually performs the difference equation calculations 
is written in 8086/8087 assembly language. Admittedly, assembly language programming is a somewhat tedious 
task; unfortunately, high-level languages (such as Pascal or Fortran) generate machine-level programs that are 
much longer, and hence much slower to execute, than an equivalent assembly language program. Given the com- 
putational constraints imposed by the PC, assembly language implementation of the digital PD algorithm is the 
only reasonable approach. The maximum sampling rate achieved for this software is 2.86 kHz (T = 350 pS). 


SECTION 5 - Frequency Dependence of Stiffness and Damping Characteristics 

In Section 3 it was shown that implementation of ideal PD control produced a system in which the stiffness 
and damping properties of the magnetic bearing could be determined by choice of two controller parameters K p 
and K d . In practice, however, such an ideal controller cannot be built. Thus, we investigate the effects of finite- 
bandwidth controllers on the bearing’s stiffness and damping properties [Af^(w) and C eq (co)]. 

Figures 8A, 8B, 9A and 9B present calculated values of K eq (wi) and C eq ( <D) for the analog controller as func- 
tions of the radial excitation frequency encountered at the bearing location. The ideal curves represent the 
expected behavior of the ideal PD controller postulated in equation (13); the actual curves show the calculated 
stiffness and damping properties of the controller actually tested for this work [equation (16)]. Note that at low 
frequencies, the ideal and actual controllers behave quite similarly. As the excitation frequency is increased, 
however, the high frequency limitations of the actual controller are no longer negligible, and their presence starts 
to degrade the performance of the magnetic bearing. In addition, the independent behavior of the ideal stiffness 
and damping properties as functions of K p and K d is not fully realized with the actual controller: varying K g has a 
slight effect on the damping curve (Figure 8B), and changes in K r can be seen in the stiffness properties (Figure 
9A). 

Similar stiffness and damping curves are presented in Figures 10A, 10B, 11A and 11B for the digital PD 
controller, although for these curves, definition of an ‘ideal’ controller is a somewhat arbitrary task. In general, 
higher sampling rates result in improved controller performance, so our ideal digital controller is defined as one 
which has a sampling rate of 1 MHz (roughly 400 times faster than the actual digital controller presented here). 
Once again, the bandwidth limitations of the actual controller are clear: bearing performance is degraded at high 


57 



frequencies, and the relationships between the controller parameters ( K p , K d ) and the stiffness and damping pro- 
perties are not completely independent. 

While there is clearly some deviation from ideal controller behavior at high frequencies, we note that both 
the analog and digital controllers behave in a very similar fashion, in spite of the fact that the resistor-capacitor 
network that makes up the analog control system is obviously quite different from the regularly re-evaluated 
difference equation in the digital control system. This similarity results from the imperfect high frequency 
derivative calculations performed by both controllers: analog differentiator performance is limited by the pres- 
ence of the high frequency poles mentioned earlier, and the straight-line derivative approximation used in the 
digital control system becomes progressively less accurate as the input signal frequency is increased. Note that a 
faster sampling rate would improve the 'high frequency performance of the digital controller, as increasing the 
roll-off points of the poles associated with the analog differentiator would increase the usable bandwidth of the 
analog controller. 

SECTION 6 - Experimental Results 

The test rig used to evaluate bearing performance is shown in Figure 12. Three shaft masses were mounted 
on the 1/2 inch diameter flexible shaft as shown; the magnetic bearing journals each weigh 5,12 N (1.15 lb), and 
the center mass weighs 7.88 N (1.77 lb). For all of the experimental results gathered, the inboard bearing was 
controlled by the analog controller with parameter settings ^=2.0 and K r = 4.0. Two sets of tests were performed: 
the stiffness of the bearing was varied while the damping was held constant, and the damping was adjusted while 
the stiffness remained unchanged. These tests were run for both the analog and digital PD controllers, and the 
results presented represent the vertical vibration response at the outboard bearing (outboard horizontal responses 
were very similar to the vertical results). 

The effect of varying the proportional feedback (stiffness) for both controllers is presented in Figures 13A 
and 13B. The curves represent the vertical amplitude response of the rotor as observed at the outboard bearing 
(similar results were obtained for the inboard bearing and are not presented here). The curves clearly indicate that 
bearing stiffness is increasing with greater amounts of proportional feedback: for frequencies below the first criti- 
cal speed (about 1800 RPM), the amplitude of vibration is reduced as the proportional feedback is increased. 

At higher frequencies, there are two other peaks (critical speeds) at about 3700 RPM and 5500 RPM. 
Changes in the bearing stiffness have little effect on the second critical because the rotor has a conical mode 
shape for this critical. It is nearly a rigid body mode shape with the ends rotating 180 degrees out of phase, and 
the center is a node point, so the bearing damping effectively keeps this peak small. The third critical is quite dif- 
ferent. It is a bending mode with the end masses (bearing magnetic laminated disks) out of phase with the center 
mass. Apparently, increasing the bearing stiffness increases the peak vibration level, a result which is counter to 
intuition. What seems to be happening is that the magnetic force is applied at the disk location and probably does 
decrease the shaft motion at the bearing location, but the sensor is located a small distance away from the disk 
center and reports an apparent increase in the peak vibration level. Thus, the measured amplitude of motion 
increases with increasing K p . 

Rotor responses in the presence of various amounts of bearing damping are presented in Figures 14A and 
14B. For both the analog and digital controllers, the measured results confirm the expectations for this test: as 
the bearing damping is increased by increasing the derivative feedback, the peak amplitude responses at the criti- 
cal speeds are significantly reduced, and the critical frequencies themselves are moved to higher frequencies (this 
last effect is most noticable at the third critical speed, which varies from 4600 RPM to 5600 RPM as K d is 
increased from 2.0 to 3.5 via the digital controller). 

SECTION 7 - Conclusions 

A PC-based digital control system has been designed and tested for controlling a radial magnetic journal 
bearing. Implementation of a proportional-derivative control algorithm on the digital computer has been shown 
to yield similar results to those obtained with an analog PD controller. For both controllers, varying the propor- 
tional gain coefficient altered the stiffness of the magnetic bearing, and varying the derivative gain changed the 
effective damping of the magnetic bearing. 
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The frequency dependent nature of the effective stiffness and damping properties of the magnetic bearing 
has been investigated. High frequency poles present in the actual analog controller introduced phase-lag into the 
control signal, and degraded operation of the magnetic bearing system at higher frequencies. Computational 
delays and zero-order-hold effects in the digital controller were seen to have effects similar to those produced by 
the high frequency poles present in the analog system. 

Both the analog and digital controllers were effective at altering the dynamic properties of the bearing under 
test. Judicious choice of the controller parameters (K g & K r for the analog controller, K p & Kj for the digital 
controller) allows significant alteration of the vibration amplitudes and critical speeds of the rotor system. 

The introduction of digital control into a magnetic bearing system offers the potential for evaluation and 
testing of algorithms much more complex than would be possible with an analog controller. For example, addi- 
tion of an integrator to the controller creates a system which rejects static loads on the rotor. This addition is a 
somewhat tedious chore with an analog system, but with a digital controller it can be accomplished by merely ins- 
talling a different software program on the computer. Also, complex digital control algorithms which can ‘shape’ 
the K eq ( ©) and C eq (i o) curve to meet certain specifications can be developed and tested. Analog realization of 
such algorithms would be very difficult, if it were even possible. Thus, digital control provides bearing designers 
a flexible tool with which they can design bearings exhibiting dynamic properties unlike those of any passive 
bearings. 
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APPENDIX A - Parameters for Magnetic Bearing Control System 
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Figure 1 - Magnetic Bearing Configuration 
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Figure 2 - Forces on Magnetic Bearing Journal 













Figure 5 - Single Degree-of-Freedom Bearing Model 


GAP 

VOLTAGE 



OUTPUT 

LOW-PASS V ^ 

FILTER VOLTAGE -current 

AMPLIFIER 


Figure 6 - Analog Controller (Block Diagram) 
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Figure 12 - Experimental Apparatus 
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Figure 13B - Bearing Vertical Displacement - Varying Proportional Feedback 
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ROTOR SPEED (RPMX10 3 ) 


Figure 14A - Bearing Vertical Displacement - Varying Derivative Feedback 
(Analog Controller) 
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Figure 14B - Bearing Vertical Displacement - Varying Derivative Feedback 
(Digital Controller) 






